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ABSTRACT 



Or 



\q ' Context. IRAS observations show the existence of a correlation between the infrared luminosity L IR and dust temperature T A in star-forming 
) , galaxies, in which larger L IR leads to higher dust temperature. The Lir-^j relation is commonly seen as reflecting the increase in dust temper- 
ature in galaxies with higher star formation rate (SFR). Even though the correlation shows a significant amount of dispersion, a unique relation 
has been commonly used to construct spectral energy distributions (SEDs) of galaxies in distant universe studies, such as source number count- 
ing or photometric redshift determination. 
O ' Aims. In this work, we introduce a new parameter, namely the size of the star-forming region r IR and lay out the empirical and modelled relation 

, between the global parameters Lir , Ta and r IR of IR-bright non-AGN galaxies. 
(/} ■ Methods. IRAS 60-to-100/jm color is used as a proxy for the dust temperature and the 1.4 GHz radio contiuum (RC) emission for the infrared 
spatial distribution. The analysis has been carried out on two samples. The first one is made of the galaxies from the 60 fim flux-limited IRAS 
^ , Revised Bright Galaxy Samples (RBGS) which have a reliable RC size estimate from the VLA follow-ups of the IRAS Bright Galaxy Samples. 
K~Ji ' T ne secon d is made of the sources from the 170 /im ISOPHOT Serendipity Sky Survey (ISOSSS) which are resolved by the NRAO VLA Sky 
K> Survey (NVSS) or by the Faint Images of the Radio Sky at Twenty-cm survey (FIRST). 
^ Results. We show that the dispersion in the L^-T^ diagram can be reduced to a relation between the infrared surface brightness and the dust 
■ - - ' temperature, a relation that spans 5 orders of magnitude in surface brightness. 

Conclusions. We explored the physical processes giving rise to the Eir-Tj relation, and show that it can be derived from the Schmidt law, 
which relates the star formation rate to the gas surface density. 

Key words. Galaxies: fundamental parameters - Galaxies: starburst - Infrared: galaxies - Radio continuum: galaxies 



1. Introduction 

The sky survey by the IRAS satellite (Neugebauer et al. 1984 1 
led to the discovery of strong connections between global pa- 
rameters of galaxies in the local universe. Among them, the 
60-to-100yum flux density ratio /?(60/100) versus L IR IRAS dia- 
gram (Soifer et al. 1987 1 exhibits a large dispersion. The quan- 
tities L IR and /?(60/100) of quiescent and starburst galaxies are 
fundamental parameters for the study of star formation. The 
first one represents the energy absorbed and reprocessed by 
dust and is related to the star formation rate. The second pa- 
rameter traces the dust temperature 7a and also provides an 
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estimate of the star formation efficiency as defined by the SFR 
per unit of gas mass (Young et al. 1986 Chini et al. 1992 1. 

The first attempts to relate the infrared surface brightness 
(which we also refer to as infrared compactness) to the dust 
temperature were hampered by the lack of sufficient infrared 
spatial resolution. Devereux ( 1987 1 used ground-based small- 
beam observations at 10 //m and compared them to the large- 
beam 12yum IRAS flux densities to estimate the compactness of 
optically bright galaxies. He showed that the ratio between the 
small-beam and large-beam flux densities is correlated with the 
global IRAS 12-to-25/vm flux density ratio which also traces 
the dust temperature. The main limitation of this work is the use 
of a rough compactness estimator that can not be easily related 
to physical parameters and thus no concluding relationship was 
derived. 
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Other similar attempts showed that the optical surface 
brightness (Phillips & Disney 1988 1 and the infrared sur- 
face brightness derived from H„ effective area (Lehnert & 
Heckman 1996 1 of IR bright galaxies increase with /?(60/100). 
However, the quantitative understanding of these correlations 
is not straightforward because the optical surface brightness re- 
sults from stars that may not be related to the dust emission and 
because in the second study, in addition to the small size of the 
sample (32 galaxies), the authors used H ff maps to estimate the 
star-forming region size without applying an extinction correc- 
tion which turns out to be crucial (Chanial et al., in prep). 

Wang & Helou (1992i made use of the extinction-free 
1 .4 GHz radio continuum (RC) size estimators to show that the 
infrared luminosity is not proportional to the galaxy physical 
area. They also showed an empirical relation between the mean 
RC surface brightness Xrc and luminosity Lrc, but they did not 
consider the IR color variations within their sample which, as 
we will show in this article, is directly related to the scatter in 
the Xrc-Lrc relation. 

More recently, Roussel et al. (2001 1 studied a sample 
of galaxies mapped by the ISOCAM camera on board ISO 
(Cesarsky et al. 1996 1 and found a correlation between the 
15-to-7yum flux density ratio and the 15/iin effective surface 
brightness. Their sample only consists of quiescent spirals of 
moderate infrared luminosities and the authors only considered 
the circumnuclear region. 

In this paper, we extend these approaches by (1) studying 
statistically larger samples, (2) applying a flux-limited selec- 
tion criterion at 60 fim, (3) considering fundamental parameters 
representative of the whole galaxy, and (4) using an extinction- 
independent estimator of the star-forming region size, the radio 
continuum angular size. 

Section 2 describes the global parameters used in this pa- 
per and Sect. 3 presents the galaxy samples that are used in 
Sects. 4 and 5 for the analysis of the luminosity-temperature 
and compactness-temperature relations. The latter is modelled 
in Sect. 6 and discussed in Sect. 7. 



2. Parameter definitions and estimations 

2.1. The dust temperature 

The temperature of the dust in a given galaxy is subject to spa- 
tial variations (Dale et al. [T999 1 that show a decrease from the 
center outwards. So, the attempt to describe the global ther- 
mal emission of the dust with one or two modified blackbod- 
ies may be questioned. However, submillimeter observations at 
450/rni and 85 //m ( Dunne et al. l2000l Dunne & Eales l200T1 
Vlahakis et al. 2005) have shown that the global emission at 
these wavelengths, which trace cold dust, does correlate with 
the global T?(60/100), which traces warmer dust. This finding 
suggests that inner and outer emissions may be bound together. 
Such a property would likely be the manifestation of the star 
formation regulation occurring on a global scale, as it appears 
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in the Schmidt law (I1959t or in relations which involve the 
galaxy rotation curve such as the Toomre's stability criterion 



: 19641 and the law by Kennicutt ( 1998 1. 



30 35 40 45 
T d (60/100) [K] 



Fig. 1. Temperature of the modified blackbodies (jS=1.3) fit- 
ted to the 60 and 100 fim IRAS flux densities (abscissae) and 
to additional far-infrared and submillimeter observations (or- 
dinates). The SLUGS 450 fim subsample and the SINGS sub- 
sample are represented with filled and open circles. The dashed 
line is the unity line and the average uncertainties are shown in 
the upper left corner. 



The often used two-component model does not reflect the 
fact that a single free parameter such as /?(60/100) suffices to 
describe the SEDs of galaxies whose infrared emission arises 
from star formation (Dale et al. 2001] Dale & Helou l2"0"0"2l . 
As a matter of fact, Serjeant & Harrison (2005) constructed 
a library of two-components IR templates, but to do so, they 
parametrized the temperatures and the relative weights with 
T?(60/100). By contrast, using a single blackbody with a con- 
stant emissivity index certainly does not provide the most accu- 
rate description of the SED, but it may supply a more profound 
insight on theglobal state of the dust in a galaxy. 

We thus define for this paper the effective dust temperature 
T& of a galaxy as the temperature of the modified blackbody of 
fixed emissivity index /3 which best fits the galaxy rest-frame 
SED between 50 and 1000 /mi. The SCUBA Local Universe 
Galaxy Survey (Dunne et al. 2000 1 is to date the largest ho- 
mogeneous set of infrared bright galaxies observed at 850 fim. 
The authors derived a mean emissivity index <y3>=1.3 with a 
standard deviation of 0.2. A subsample was then observed at 
450 fim (Dunne & Eales 2001 ) but the authors did not publish 
their single-temperature analysis, althought they do state that 
this model cannot be ruled out, and instead favoured the anal- 
ysis of the two-component model which offers one additional 
free parameter. By excluding galaxies with an active galactic 
nucleus (AGN) from their 450/mi sample and by taking the 
flux densities tabulated in their article, we found a similar value 
of <y8>=1.38 with a sample standard deviation of 0.17. Since 
only 2 out of 25 galaxies have a reduced ^>2, the single- 
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Name 


Morphology 


D 


logL, R / v (60^m) /.(lOOyum) T d (60/100) 


T A 


(1) 


(2) 


(3) 


(4) (5) (6) (7) 


(8) 



UGC 903 


Sc 


33 


.16 


10 


.32 


7 


.78 


15 


.45 


35.5 d 


= 2 


.7 


33.5 = 


z 1 


.0 


NGC 958 


SBc 


76 


.36 


11. 


.00 


5 


.85 


15 


.08 


31.5 d 


2 


.2 


30.1 = 


z 0. 


.7 


UGC 2369 S 


Sbc 


121 


.90 


11. 


.39 


8 


.07 


11 


.18 


40.4 d 


z 3 


.3 


40.7 = 


z 1. 


.5 


UGC 2403 


SBa 


53 


.78 


10. 


.68 


7 


.72 


12 


.06 


38.3 d 


z 3 


.1 


36.2 = 


z 1. 


.2 


NGC 2856 


Sc 


40 


.08 


10. 


,32 


5 


.73 


10 


.15 


37.0 = 


z 2 


.9 


36.2 = 


z 1 


.3 


NGC 2990 


Sc 


47 


.18 


10. 


.43 


5 


.16 


9 


.61 


35.5 d 


z 2 


.7 


34.2 = 


z 1 


.1 


UGC 5376 


Sd 


31 


.28 


10. 


.10 


5 


.36 


10 


.41 


34.7 d 


z 2 


.6 


33.1 = 


z 1. 


.0 


Arp 148 E 


Irr 


143 


.20 


11. 


.45 


6 


.38 


10 


.30 


38.9 d 


z 3 


.1 


38.2 = 


z 1. 


.4 


Zw 247.020 


Sa 


107 


.70 


11. 


.16 


6 


.01 


8 


.47 


40.8 d 


z 3 


.4 


43.9 = 


z 2. 


.0 


I Zw 107 


Pair 


168 


.70 


11. 


.70 


9 


.02 


10 


.00 


46.4 d 


z 4 


.2 


44.6 = 


z 1 


.9 


NGC 5962 


Sc 


32 


.17 


10. 


.41 


8 


.93 


21 


.82 


32.3 d 


= 2 


.3 


31.3 = 


z 0. 


.8 


NGC 6052 


Pair 


70 


.42 


10. 


.86 


6 


.79 


10 


.57 


38.3 = 


= 3 


.1 


36.8 = 


= 1. 


.3 


NGC 6181 


SABc 


30 


.70 


10. 


.36 


8 


.94 


20 


.83 


32.7 d 


z 2. 


.4 


33.3 = 


= 1. 


.0 


NGC 7541 


SBbc 


30 


.09 


10. 


66 


20. 


08 


41 


.87 


34.5 d 


z 2. 


.6 


34.4 = 


= 1. 


.0 


NGC 520 S 


Sa 


30 


.22 


10. 


.78 


31. 


.52 


47. 


.37 


39.1 d 


z 3 


.2 


40.7 = 


= 1. 


.3 


UGC 2982 


SABc 


67 


.57 


10. 


.97 


8 


.39 


16 


.82 


34.7 d 


z 2 


.6 


35.4 z 


= 1 


.0 


NGC 2623 


Sa 


77 


.43 


11. 


.44 


23. 


.74 


25 


.88 


46.2 d 


z 4 


.2 


49 .1 z 


= 1 


.8 


NGC 3110 


SBb 


73 


.48 


11. 


.17 


11. 


.28 


22 


.27 


34.7 d 


z 2 


.6 


36.7 z 


= 1. 


.1 


IRAS 10173+0828 




198 


.70 


11. 


.63 


5. 


.61 


5 


.86 


49.5 d 


z 4 


.6 


46.9 d 


: I, 


.7 


IRAS 10565+2448 




176 


.30 


11. 


.87 


12. 


10 


15 


.01 


43.2 d 


z 3 


.7 


47.2 d 


z I. 


.6 


IRAS 12112+0305 NE 




292. 


.50 


12. 


.14 


8 


.18 


9 


.46 


48.1 d 


z 4 


.3 


48.0 d 


: I. 


.8 


UGC 8387 


Irr 


99 


.99 


11. 


.55 


17. 


.04 


24 


.38 


36.2 d 


z 2 


.8 


41.3 d 


: I. 


.2 


Zw 049.057 


Irr 


59 


.06 


11. 


.20 


21. 


.89 


31 


.53 


40.0 d 


z 3 


.3 


40.9 d 


: I. 


.3 


NGC 7592 


Pair 


95 


.13 


11. 


.17 


8 


.05 


10 


.58 


42.0 d 


z 3 


.5 


38.7 d 


: I. 


.5 


NGC 7714 


Sb 


38 


.16 


10. 


.50 


11. 


16 


12 


.26 


44.8 d 


z 4 


.0 


43.2 d 


z I. 


.8 



Columns: (3) distance in Mpc; (4) infrared luminosity in L© derived from the 60 and 100 [im IRAS flux densities by assuming a modified blackbody of emissivity index /3=1.3; (6) & 
(7) IRAS flux densities in Jy; (7) effective dust temperature and lcr uncertainties in K, 1 .3; (8) effective dust temperature and lcr uncertainties in K derived from the additional data 
points tabulated in Dunne & Eales i 2001 '. 



temperature model provides a reasonable fit in most cases and 
we adopted the fiducial value /?= 1 .3 in this paper. 



to validate the use of /?(60/100) to estimate the effective dust 
temperature in our study. 



The next step is to check that it is possible to easily asso- 
ciate an estimator to the previously defined effective dust tem- 
perature, the estimator of choice being the color /?(60/100) due 
to its matchless availability in the local universe. We have al- 
ready noted that the T?(60/100) parameter can characterize to 
some extent the whole infrared SED of normal and starburst 
galaxies, but its precision still has to be determined. For two 
sets of non-AGN galaxies described in Tables ^ ar, d ID we 
compared the effective dust temperatures exclusively derived 
from the 60 and 100 fim IRAS flux densities to the temperatures 
that are derived from additional far-infrared (FIR) and sub- 
millimeter observations. The two sets are the SCUBA SLUGS 
450 /im subsample (Dunne & Eales 2001 1 and the galaxies ob- 
served by Spitzer at 70 and 160 /jm from the SINGS catalog 
(Kennicutt et al. l2"0"0"3l Dale et al. l2005> for which submillime- 
ter observations (800 or 850/vm) were available in the litera- 
ture. We completed the second set with far-infrared flux densi- 
ties mainly from ISOPHOT (Lemke et al. 1996 1 and the Kuiper 
Airborne Observatory. The results are shown in Fig.^and con- 
firm /?(60/100) as an estimator of the effective dust tempera- 
ture. The standard deviation is 2.8 K, which is small enough 



2.2. The IR luminosity 



The infrared luminosity is estimated from the bolometric lumi- 
nosity of the best-fitting blackbody modified by a A. - 3 emis- 
sivity function. As a consequence, it is a measure of the IR 
emission of the big dust grains in thermal equilibrium and 
it excludes the emission from smaller dust grains stochasti- 
cally heated such as the polycyclic aromatic hydrocarbons. 
Distances for Hq = 75kms~ 1 are taken from, in decreas- 
ing order of precedence: the Catalog of Neighboring Galaxies 



(Karachentsev et al. 2004 1, the Revised Bright Galaxy Sample 
(Sanders et al. 2003 1, the Surface Brightness Fluctuation survey 
(Tonry et al. l200Tl . the Nearby Galaxies Catalog (Tully[l988l 
and the NASA/IPAC Extragalactic Database (NED) whose red- 
shifts were corrected for the Virgo inflow. IRAS flux densi- 
ties are taken from, by order of precedence: the RBGS, the 
Large Optical Galaxies Catalog (Rice et al. 1988 1, the Point 
Source Catalogue with redshift (PSCz, Saunders et al. 2000 1, 
the Point Source Catalog (PSC, Beichman et al. U988> and the 



Faint Source Catalog (FSC, Moshir et al. [19921 
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Table 2. /?(60/100) as an effective temperature estimator: the Spitzer SINGS subsample. 



Name Morphology D \o\ 


gL m / v (60//m)/ v (70/im) / v (100/xm) / v (160/zm) A 


jM) 


Refs. 


r d (60/ioo) 




1 1 1 


n\ 
\A) 






(j) 






{"> 


\ J ) 


(10) 


(11) 


(12) 


(13) 


NGC 337 


SBcd 


21.6 


10.04 


9.07 


8.83 


20.11 


18.30 


850 


U.JJ 


1 


39 £ -u 
jZ.O ± 


9 A 




-1- 1 1 
±1.1 


NGC 2798 


SBa 


27.8 


10.52 


20.60 


14.70 


29.69 


18.45 


850 


0.19 


I 


3Q A -A- 


3 9 

J.L 


jo.j 


-1- 1 fx 

± 1 .o 


NGC 2976 


Sc 


3.6 


8.68 


13.09 


16.99 


33.43 


46.81 


850 




j 


jU.o ± 


L.L 


Ly.V 


-1- 1 1 
± L.L 


NGC 3190 


Sa 


24.1 


9.80 


3.33 


4.34 


9.84 


13.19 


850 


n iq 

U.17 


1 


9Q 3 -u 
Ly.j ± 


z.u 


Lo.j 


_i_ n q 


Mrk33 


Sm 


26.8 


9.83 


4.79 


3.34 


5.49 


3.46 


850, 850 


05 04 


2. 1 


A!X Q -4- 


3 Q 
J.y 


J7.U 


-i- 1 ^ 

± L .J 


NGC 3521 


SABb 


6.8 


9.81 


49.19 


49.85 


121.80 


206.70 


850 


2.1 1 


1 


jl.j ± 


9 9 
L.L 


Zo.o 


-A- 1 9 
± L.L 


NGC 4254 


Sc 


15.3 


10.39 


37.46 


39.02 


91.86 


131.80 


160, 350, 360 


78 7.8 16 


J, T - , J 


j 1 .0 ± 


L.L 


jj.U 


± L.VJ 


















450, 800, 850 


3 8 fi 101 


4, 4, 1 










NGC 4321 


SABb 


15.2 


10.25 


26.00 


32.28 


68.37 


128.40 


160, 850 


46 0.88 


3, 1 


30.7 ± 


2.1 


30.5 


± 1.2 


NGC 4536 


SABb 


14.9 


10.15 


30.26 


22.49 


44.51 


54.39 


850 


0.42 


] 


39.1 ± 


3.2 


35.1 


± 1.6 


NGC 4631 


SBcd 


7.7 


10.08 


85.40 


98.78 


160.10 


269.00 


180,450,450 


121,25.06, 30.7 


5,6, 1 


35.0 ± 


2.7 


30.7 


±0.7 


















450, 850, 850 


18,5.253,5.73 


5,6,1 


























850, 870 


1.89,3.78 


5,7 










NGC 5195 


SBO 


7.7 


9.35 


15.22 


10.85 


31.33 


12.34 


850 


0.26 


1 


33.7 + 


2.5 


35.7 


± 1.5 


NGC 5713 


SABb 


26.7 


10.54 


22.10 


17.23 


37.28 


34.77 


180,450, 800 


16, 0.889,0.102 


5,8,8 


36.8 ± 


2.9 


38.4 


± 1.2 


















850, 850, 850 


0.359,0.57,0.43 


9,1,5 










NGC 5866 


SO-a 


15.3 


9.65 


5.26 


6.66 


16.98 


16.53 


180,450, 850 


10.5,0.79,0.14 


5,1,1 


28.3 ± 


1.9 


31.9 


±0.8 


NGC 6946 


SABc 


5.9 


10.08 


129.80 


177.90 


290.70 


498.40 


60, 60, 100 


165, 115.5,338 


10, 11, 10 


32.5 ± 


2.4 


31.7 


± 0.7 


















160, 200, 200 


450, 330, 743 


12, 12, 10 


























200,450, 850 


365.8,18.53,2.98 


11, 1, 1 










NGC 7331 


Sb 


13.1 


10.34 


45.00 


56.49 


110.20 


164.10 


60, 100, 200 


42.9, 120,243 


10, 10, 10 


31.4 ± 


2.2 


28.9 


±0.8 


















450, 850 


20.56, 2.11 


1, 1 










NGC 7552 


Sab 


21.4 


10.86 


77.37 


45.40 


102.90 


86.65 


850 


0.8 


1 


40.9 ± 


3.4 


36.4 


± 1.6 



Columns: (3) distance in Mpc; (4) infrared luminosity in L© as described in TablefTl (5) & (7) IRAS flux densities in Jy; (6) & (8) Spitzer flux densities in Jy; (9) & (10) Flux densities 
in Jy and references for the additional far-infrared and submillimeter data: [1] Dale et al. 120051 . [2] Hunt et al. 120051 . [3] Stark et al. 119891 . [4] Eales et al. 119891 , [5] Bendo et al. 
I^OT^l . [6] Stevens et al. 120051 . [7] Dumke et al. [8] Chini et al. [19951 . [9] Dunne et al. I2TOI [10] Alton et al. tTffSgl . [11] Tuffs & Gabriel EuTBI . [12] Engargiola tTWTl : 

(12) effective temperature and ltr uncertainties in K derived from the 60 and 100 /im IRAS flux densities; (13) effective temperature and lcr uncertainties in K derived from all data 
points. 



2.3. The star-forming size 



2.4. Infrared surface brightness 



A well-known tracer of star formation, for which a wealth of 
large and high-resolution catalogs are available, and that is un- 
affected by dust extinction is the radio continumm. Globally, 
the luminosity of non-AGN galaxies at 1 .4 GHz correlates with 
the FIR luminosity over 4 orders of magnitudes (Condon 1992 
Yun et al. 2001 1. Spatially, Chanial et al. (in prep.) show that 
the sizes of non-AGN galaxies inferred from radio continuum 
observations correlate tightly with the size inferred from FIR 
and CO maps, even for infrared luminous galaxies, in the CO 
case. 

In this paper, the star-forming size is estimated (Chanial et 
al. in prep) as 



(0.86 ± 0.05) r RC , 



(1) 



where r RC is the HWHM of the deconvolved RC emission (see 
Sect. 3). The HWHM is chosen to be along the major axis, to 
account for the galaxy inclination at a first order. 



We define the observed infrared surface brightness by the for- 
mula 



Sir — 



(2) 



The factor i has been introduced to provide a better esti- 
mate of the IR emission from the effective area nr^ R , where r m 
is the observed HWHM maximum of the emission along the 
major-axis, because for axisymmetric gaussian profiles (face- 
on galaxies case), the half-light radius r(|)= V21n 2cr is equal 
to the HWHM. 



3. Sample definitions 



Soifer & Neugebauer ( 1991 ) showed that the complete lOOjt/m 
flux-limited subsample of a complete 60 fim flux-limited sam- 
ple has a colder average dust temperature, which implies that 
the Li R -Ta relation is biased by the wavelength at which a sam- 
ple is selected (see also Blain et al. 2004). Furthermore, cold ul- 
traluminous infrared galaxies falling off this relation have been 
discovered in the 170yum FIRBACK survey (Chapman et al. 
2002) and in submillimeter surveys (Chapman et al. 2005 1. 
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To test whether the infrared compactness-temperature rela- 
tion is subject to such a bias, we performed our analysis on two 
samples of non-AGN galaxies selected at 60 and 170/mi. 

3.1. The 60 /urn selected sample 

This sample has been obtained by Chanial et al. (in prep.) 
by matching the RBGS, which is complete for extragalactic 
sources with / v (60//m)>5.24Jy and Galactic latitudes \b\>5°, 
with the two RC follow-ups around 1.4 GHz (Condon et al. 
rH%Ul IT9961 of the 60 /mi flux-limited IRAS Bright Galaxy 
Samples (Soifer et al. [T^87|[T989l and Sanders et al. fl995t . For 
each IRAS source, one or more maps have been observed with 
an angular resolution ranging from 1'.'5 to 60". Deconvolved 
major-axis FWHM are provided by both datasets and are de- 
rived from two-dimensional gaussian fits to the maps. 

Starting from this initial list, we applied several editing 
steps to ensure the reliability of the size estimates (cf. Wang 
& Helou fTffHl Meurer et al. 119971 . 

1. We retained only spatially resolved sources with angular 
sizes no lesser than half the radio beam FWHM. 

2. Because at higher angular resolutions, extended emission 
may be missed, we retained only the radio sources that con- 
tribute to more than 2/3 of the total radio continuum flux, ensur- 
ing that the angular size is representative of the whole galaxy. 
This step also ensures that no more than one radio counterpart 
is retained. 

3. For most of the IRAS sources, no more than one RC map goes 
through the two steps above and thus no more than one angular 
size estimate was deemed reliable. For the other sources, we 
adopted the mean value of the angular sizes. 

4. Contamination by an AGN has been dealt with very conser- 
vatively by excluding galaxies satistying any of the following 
optical, infrared and radio criteria: (i) position closer than 30" 
to an AGN galaxy (including LINERs) from the extensive cata- 
logue by Veron-Cetty & Veron ( 2006 1 or classification in NED 
as AGN, (ii) / v (25//m)// v (60/<m)>0.2 (de Grijp et al.fT985t 
and (iii) low IR-to-radio ratio q<1.94 (see Yun et al. 2001 1. 

The sample resulting from this selection process contains 
232 sources. They are mostly IR-dominated spirals or inter- 
acting galaxies and their global properties are summarized in 
Fig.0 

3.2. The 1 70 /urn selected sample 

The second sample has been extracted from the ISOPHOT 
Serendipity Sky Survey (Stickel et al I2004> . which is made 
of 1927 sources detected at 170 /mi by the C200 ISOPHOT 
detector during the slews between pointed observations. The 
ISOSSS covers 15% of the sky. A systematic inspection of the 
DSS (Digitized Sky Survey) snapshots of the sources and ex- 
amining the FSC and the PSC lead us to discard 52 star iden- 
tifications and 20 possible star associations, one likely plane- 
tary nebula, one extragalactic H n region, one radio source not 
associated with the ISOPHOT source and 7 entries that may 
be contaminated by a nearby radio source. We also excluded 
two offcenter duplicate entries of NGC 7331 and a source con- 



taminated by cirrus. Before cross-correlating the ISOSSS cat- 
alog with radio catalogs, we made sure that the optical candi- 
dates taken from the HYPERLEDA database within 2' of the 
ISOSSS position had an astrometry with one arcsecond accu- 
racy. This work resulted in the removal of 38 duplicates, the 
addition of 241 new sources (including 160 stars), and 515 new 
positions, 426 of which were taken from the 2MASS Extended 
Source Catalog (Jarrett et al. 2000 1 the remaining ones were 
determined on IR POSS-2 plates, which absolute astrometry 
were corrected using the USNO-B1 and GSC2.2 star catalogs. 

Then, we cross-correlated the optical candidates with 
two complementary 1.4 GHz RC surveys: the NVSS catalog 
(Condon et al. [T9981 and the FIRST catalog (April 2003 re- 
lease, Becker et al. I1995> . The former has a 45" resolution 
beam and covers 82% of the sky while the latter has a higher 
resolution (5") but a more limited sky coverage (22%). Finally, 
after cross-identification with the NED database, AGN can- 
didates were removed by applying the same methodology as 
adopted for the 60 /mi sample. At this point, the sample con- 
tains 899 remaining sources. To avoid undersampling effects, 
we only retained sources whose NVSS counterpart has an an- 
gular major axis 22'.'5<6> RC <120" or whose FIRST counterpart 
has an angular major axis 2'.'5<# RC <15". NVSS and FIRST 
major axis 8 is derived from a two-dimensional gaussian fit, 
so that sources with a complex geometry such as spirals with 
dominant H n regions or closely interacting pairs are unlikely 
to be fitted. We attempted to exclude them by rejecting those 
with a radio-optical position offset greater than \& RC - We also 
discarded the FIRST sources that did not make up 75% of the 
total radio flux density, assumed to be the NVSS one. 

In the final tally, the 170 /mi sample is made of 198 galaxies 
with reliable angular size information, of which 1 34 are taken 
from the NVSS catalog and the other 64 from the FIRST cat- 
alog. Like the 60 fim sample, the two subsamples are mostly 
made of IR-dominated spirals and interacting galaxies. Their 
global properties are listed in the online Tables 4 and 5 and 
summarized in Fig. [2] 

4. Luminosity-temperature relation 

The response of the large dust grains to the heating radia- 
tion field has been studied through the Lir-T?(60/100) diagram. 
For 60 //m flux-limited samples, it has been shown that the 
dust gets warmer as the FIR luminosity increases (Smith et al. 
IT9871 Rowan-Robinson et al. lT987l Soifer et al. ll9§9l Soifer & 
Neugebauer 1991 1. This trend is also followed by the 170/mi 
sample, but more weakly. 

It has been noted that the Lvr-/?(60/ 100) relation is affected 
by a large intrinsic dispersion, as shown in Fig. [3] for our 60 
and 170/mi samples, which is attributed to the scatter of the 
gas content (Soifer et al. 1 19871 Sanders et al. U991> . Chapman 
et al. (2003 1 carried out a phenomenological study of this dis- 
persion in the IRAS 1.2 Jy survey and showed that it is key to 
understanding the populations of flux-limited surveys. 

The systematic effect of the sample selection on the Lir- 
T?(60/100) relation, which prompted us to analyse two sam- 
ples selected at different wavelengths (Sect. [3}, is confirmed 
in Fig. |4] in which the two samples are compared in 6 IR lu- 




Fig. 2. Distribution of a) distance, b) IR luminosity, c) effective dust temperature, d) IR physical radius, e) morphology type, f) 
B-band luminosity with v=c/0.44 /mi, g) IR to B-band luminosity ratio and h) FWHM of the radio continuum emission along 
the major axis for the 60 /mi sample (shaded), the 170//m sample spatially resolved in NVSS (hatched) and FIRST (dotted). The 
mean distance is 46.2, 62.3 and 154 Mpc respectively, the mean L IR is 10 10 48 , 10 10 25 and 10 10 72 L , the mean T d is 35.01, 31.36 
and 34.64 K, the mean r m is 2.10, 5.10 and 2.04 kpc and the mean B-band luminosity is 10 98l) , 10 9 82 and 10 10 01 L Q . Morphology 
type and B-band luminosity are taken from the HYPERLEDA database. 



minosity bins. For each bin, the 170/mi sample is colder that 
the 60 /mi sample, the difference being most striking for Lir 
between 10 10 and 10 1L6 L o . 



5. Compactness-temperature relation 

To better understand star formation processes on a global scale, 
we introduced an additional parameter to the Lir-Lj analysis, 
the size of the star-forming region r IR derived from the FWHM 
of the RC profiles according to Eq. \l\ This new information 
is shown in the LiR-Ld diagram (Fig.[5]a) by representing the 
galaxies from the 60 /mi sample with circles of radius pro- 
portional to r IR . A systematic effect is apparent: at fixed Lir, 
T?(60/100) increases as the physical size decreases, reflecting 
the qualitative fact that the dust temperature increases as the 
average grain-star distance decreases. 

More quantitatively, we determined the plane 



+ (0.092 ± 0.007) 



logfl(60/100) = a + Mog— + clog — 

L Q pc 



(3) 



that best matches the observations in the three dimensional pa- 
rameter space. The coefficients (a,b,c) and their associated 
errors were derived by using bootstrap samples of the 60 /mi 
dataset and the resulting empirical relation is 

log £(60/ 100) = -0.663 ± 0.095 



log | ^J- (2.09 ±0.24) log(^ 



(4) 



We note that because the covariates r m and Lir scale as the 
distance of the galaxy and as its square, the value c/b--2 may 
be unduly favoured by the x 1 minimization. To check that this 
is not the case, we ran the bootstrap estimator on a fake 60 //m 
sample, for which the values of Lir and r IR are unchanged, but 
for which the values of /?(60/100) have been randomly permu- 
tated, so that the variate is uncorrelated with the covariates. In 
the resulting sample, the bootstrap estimation of c/b does not 
converge and we obtained a sample mean of <c/Z?>=-194.3 and 
a sample standard deviation of 891.2. In fact, its distribution is 
similar to a Cauchy (or Lorentz) distribution, for which mo- 
ments do not exist. We note that such a similarity for the un- 
correlated sample is not surprizing, because the distribution of 
the ratio of two independent normal distributions is precisely 
a Cauchy distribution. The clear difference between the distri- 
butions of the c/b values for the real and the fake 60 /mi sam- 
ple (Fig. |6jl is a strong indication that the global parameters 
/?(60/100), Lir and r m are indeed correlated. 

The resulting value of c/b is equal to -2 within its uncer- 
tainties, which implies that /?(60/100) scales as LiR/r 2 R or that 
for a given fixed temperature, the IR luminosity linearly scales 
as the IR area (Fig.|3]b). This finding is non-trivial and suggests 
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Fig. 5. a) The vs Lvr diagram for the 60pm sample. The radius of the circles is the IR radius, assumed to be proportional to 
the linear radio continuum HWHM along the major axis, b) IR luminosity vs the effective IR area nrL for two temperature bins. 
Galaxies with 28.5<r<]<30.5 K are plotted with open circles and the ones with 40<7d<41 K with filled circles. The solid lines 
corresponds to the empirical correlation given by Eq.[6]for the mean Ta of the two subsets. We note that the empirical relation 
does not exactly bisect the low-temperature subsample because of the RBGS surface brightness detection threshold (dashed line, 
from Wang & Helou [T992l which affect the completeness of the low-temperature subsample. 



that the disk geometry comes into play. Assuming a value of -2, 
we performed a linear regression (bisector method) 



log #(60/ 100) = -0.783 + 0.019 

+ (0.123 + 0.005) 1 1 LlR 
that can be written as 
T d = (23.5 ± 0.3) 



log I— j-21og Vpc 



-IR 



Lope 



-2 



0.052±0.0()2 



K, 



(5) 



(6) 



where T d and Sir are the effective dust temperature and the 
infrared surface brightness defined in Sect|2] 

For the resolved 170 //m sample, the relation is 



T A = (22.9 + 0.4) 



-IR 



Lope 



-2 



0.057±0.004 



K. 



(7) 



which is consistent with the equation obtained for the 60 fim 
sample. 

6. Modelling 

Two idealized scenarios were considered to interpret the phys- 
ical processes behind the empirical relation Eq.|6j the dust be- 
ing distributed in a single shell of radius equal to the observed 
star-forming radius r IR or in molecular clouds optically thin at 
far-infrared wavelengths. 



6.1. Single dust shell 

In this case, we assume that all the dust is distributed in an 
isothermal shell of radius R around a point-source starburst. 
Lehnert & Heckman ( 1996 1 have proposed that this model is in 
agreement with the EiR-iTd relation. 

The energy absorbed and emitted by a grain of size a and 
of emissivity 2 a bs(v) is 



£ a hs = 7ia 



Jo 

T 

Jo 



2absO)fy(v) dv 



4na z Q abs {v)nB v {v,T d )dv. 



(8) 
(9) 



Assuming that the emissivity QabsM - 1 in the UV-optical 
regime and gabsO) =* gabs(vo) (v/v Q f in the infrared regime, 
equations Eq.[8]and Eq.|9]yield 



E em /na 



= 1' 

Jo 



F v (v) dv 



kTA 4+l3 8abs(VQ) f + °° JC 



h ) 



4 



I 

Jo 



r 3 + fl 



1 



-d* 



(10) 
(11) 



where vo is a reference frequency,/? is the dust emissivity index, 
and x = hv/kTd in the expanded Planck function. 

Assuming that the dust shell is optically thick, the energy 
radiated by the central starburst is totally reemitted in the in- 
frared and the right-hand part of equation Ea.llOlbecomes 



_ 1 

4nR- ~ 4 S 



IR- 



(12) 
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Fig. 3. The rest-frame 60-to-100jum IRAS flux density ratio vs 
the infrared luminosity. The galaxies from the 60 fim (170/im) 
sample are represented by filled (open) circles. Are also plotted 
the sources from the 170/im sample for which no reliable size 
estimate is available. 

The right-hand part of equation Eq.^2 can a l so by rewritten 



E t Jna 



4<T/3 gabsOo) 



p > -3, 



(13) 



by using the equality (Fikhtengol'ts ll947l 



Jo e fl 



1 



djc = -rr(Jt) £(k) k>l,%ea>0, 



where T(z) = C ? 1 e ' df is the Gamma function and 
((z) = Ti7=i 1 /" z i s me Riemann's Zeta function, and by setting 



In Ic 4+ P 

^ = ^7^ r (4 + /?K(4 + /?), 



(14) 



for which crp=Q is equal to the Stefan-Boltzmann constant that 
relates the bolometric luminosity of a blackbody to its temper- 
ature. 

At the thermal equilibrium, £ab s = E em , and the dust tem- 
perature is related to the IR surface brightness by 



Ti 



4+/3 . 



4 



^CTjgQabsOo) 



(15) 



By taking /3 = 1.3 (Sect. 12. 1 1 and the standard value 
g abs (125/mi) = 7.5 x 10~ 4 (Hildebrand fT983l cf. reviews by 
Hughes et al. 1997 and Alton et al. 2004 1, the numerical rela- 
tion between and Eir in the single dust shell model is 



T d = 7.6 



L G pc 



0.19 



K. 



(16) 
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Fig. 4. Distribution of the rest-frame 60-to-100/im IRAS flux 
density ratio for 6 IR luminosity bins spanning 10 9 6 to 10 12 L Q . 
The shaded (hatched) histogram relates to the 60jum (170/vm) 
sample. The 170//m sample includes the sources for which no 
reliable size estimate is available. 



It is plotted in figure Fig.[7]and labelled as (la). It departs sig- 
nificantly from the observational data. 

To check that the difference is not due to the crudeness of 
the assumed dust model and to compare the work by Lehnert & 
Heckman ( 1996 1, we also used the more realistic dust model by 
Desert et al. ( 1990). This dust model is calibrated on solar inter- 
stellar medium abundances and includes the stochastic heating 
of PAHs. Although it assumes an isotropic radiation field, it can 
be used for the single shell geometry, as long as we input the ra- 
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Fig. 6. Distribution of the values c/b resulting from the re- 
gression fits of the bootstrap samples to the linear form 
logT?(60/100) = a + felogLiR/L + clogr IR /pc. For the shaded 
histogram, the real observed values T?(60/100), Lir, and T d of 
the 60 fj.m sample have been used whereas in the hatched his- 
togram, the R(6Q/ 100) values have been randomly permutated. 
The boostrap estimation of c/b only converges for the real ob- 
served sample, to the value -2.09±0.24, which is consistent 
with /?(60/100) being a function of L^/r 2 . 



diative energy density of the dust shell, u v = Sir /4c. The single 
shell model with the Desert et al. ( 1990 1 dust model and with 
a heating source scaling as a 05 star is plotted in Fig. Q and 
labelled as (lb). It shows a steep relation similar to equation 
Eq.^J which is not in accordance with our samples. 

As a result, the geometry itself of the single dust shell 
model is not satisfactory, unless drastic changes in the dust 
composition occur along the Sir sequence. The apparent agree- 
ment between the Lehnert & Heckman ( 1996) smaller sample 
and the single dust shell model is likely explained by the fact 
that the star-forming sizes were estimated from H a maps uncor- 
rected for dust attenuation. Because extinction preferentially 
affects high compactness regions, their maps likely missed the 
central nuclei of the most compact starbursts, leading to an 
overestimation of the star-forming sizes and an underestima- 
tion of the infrared surface brightness. 

6.2. UV-optical-thick & FIR-thin molecular clouds 

We assume that the infrared emission is from a disk of radius 
R and that the dust giving raise to this emission is isothermal. 
Several dust configurations could yield an isothermal dust pop- 
ulation such as thin shells around young star clusters or dust in 
cirrus exposed to a uniform interstellar radiation field. 

We also assume that every dust grain radiates as a black- 
body modified by a A~P emissivity function. In the case in 



which the medium is transparent in the far-infrared but opaque 
in the optical-UV, the luminosity L IR radiated by the dust is 
proportional to the total mass of dust M d and more specifically 
(HildebrandEH 



F v (v) = K(v)M d 



B v (v, Tj) 
D 2 ' 



(17) 



where F v is the flux density, k(v) is the absorption mass coeffi- 
cient assumed to be equal to k(vq) {v/v^f and D is the distance. 
The integration of this relation over the frequencies leads to 



= 4nD 2 

Jo 



F v (v)dv = 4cr/^M d r d 4 ^ 



(18) 



where crp is the "generalized" Stefan-Boltzmann constant 
which we introduced in Eg. 1141 

Setting the dust-to-gas mass ratio 77d = Md/M gas , the total 
luminosity can be related to the total gas mass M eas by 



A K ^ \* T 4+/3 
: 4 O-p—g- T] d M g:is T d P . 

^0 



(19) 



The size of the star-forming region is introduced in our 
analysis through the Schmidt law that non-linearly relates the 
star formation rate to the gas mass surface density with a power 
index that reliably departs from unity. This law has been most 
accurately determined by Kennicutt ( 1998 ) as 



£sfr 



M Q yr~' kpc 



.-2 



= 2.5 10 



,-4 



-gas 



M pc- 



1.4 



(20) 



where E gas is the gas mass surface density inside the radius R. 

Kennicutt ( 1998) derived the star formation rate from the 
total infrared luminosity 



SFR 

M G yr- 



r 



5.8 10 9 L G 



(21) 



by assuming that the dust absorbs and reprocesses all the in- 
trinsic star light (factor y=l). 

Because Eqs.^]and|^are linear, both of their sides can 
be divided by ttR 2 and expressed in terms of surface densi- 
ties and surface brightness. By substituting Esfr = SFR/nR 2 
from Eq. |^ in Eq. |20| we can derive the gas mass surface 
density as a function of the infrared surface brightness. Then, 
S gas = Mg^/nR 2 can be eliminated from Eq. 1191 and we obtain 
the following relation between T& and Eir 



T 6 



4+fi 



6.27 10" 



crp /c(vo) rjd y 



1/1.4 



L G pc 



0.4/1.4 



(22) 



Assuming that the interstellar medium is optically thick, 
y-\ and by adopting standard values of the parameters in- 
volved in our model, K(\25fim) = 1.9 m 2 kg 1 (Hildebrand 
[T983t . 7/a = 1/350 (see Sanders et al. fl99ll Bendo et al. 120031 
and/3= 1.3, Eq.l22l gives 



T d = (22.9 + 0.9) 



L G pc 



0.054±0.013 



K. 



(23) 
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Fig. 7. The effective dust temperature 7^ vs. the IR surface brightness Sir for the 60 jim sample (filled circles) and the 170/mi 
sample (open squares and triangles for the sources with reliable NVSS and FIRST angular sizes respectively). The solid thick 
line shows the empirical relation Eq.|6] The dashed line (la) represents (la) the single dust shell model with an analytical dust 
model (Sect. IrTTI Eg. 1 1 61. and the dot-dashed line (lb) the single dust shell model coupled with the Desert et al. ( 1990 > dust 
model. The long-dashed line (2) is the model involving the Schmidt law (Sect. l6~2l Ea.l23> and the dotted line (3) is the expected 
relation for a blackbody (Ea.l24>. The average uncertainties in 7a and Sir are shown based on errors of 10 per cent for/ v (60jum), 
15 per cent for / v (100/im) and 20 per cent for r IR . 



This equation is plotted in Fig.[7]and labelled as (2). It is 
in very close agreement over 5 orders of magnitude with the 
empirical relation Eq. [6] obtained from the 60 yum sample of 
galaxies and is in accord with the 170^m sample, though we 
should bear in mind that the actual position of the theoretical 



relation could be translated along the Tj-axis because of the 
uncertainties in k(vq) and r^. 
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7. Discussion 

7.1. Selection biases 

Sources with 10 10 <Lir<10 10 5 L e are on average colder in the 
170yum NVSS-resolved subsample (<r d >=30.5 K) than in the 
60//m sample (34.0 K). The difference is not the result of the 
additional selection criteria on the NVSS sources, because their 
average effective dust temperature is comparable to the whole 
170yimi sample in this luminosity bin. The temperature differ- 
ence arises from the wavelength selection that affects the L IR - 
diagram, as discussed in Sect. 0] On the contrary, the Ei R - 
Td relation is followed by both the 60 //m and 170yum selected 
samples and we conclude that the latter relation is not sensi- 
tive to the selection wavelength, and reflects more fundamental 
physical processes in play in quiescent and starburst galaxies. 

Another selection bias that could affect the £rR-7d rela- 
tion is the Malmquist bias involved in the determination of the 
Tully-Fisher relation for spiral galaxies (Tully & Fisher 1977) 
and the Fundamental Plane of elliptical galaxies (Dressier et 
al. [T9871 Djorgovski & Davis, \TWi\ . This bias would un- 
derpopulate galaxies with small physical sizes, which would 
lead to an underestimation of the effective dust temperature of 
the most compact galaxies. This effect may be present in the 
170yum samples for which the angular size is censored and to a 
lesser extent in the 60 fim sample. However, out of the 75 IRAS 
sources observed in Condon et al. ( I1990I IT996 > with an angu- 
lar resolution higher than or equal to 1'.'8, only 9 unresolved 
sources would have passed the selection criteria of the 60 fim 
sample described in Sect. 13. II 

Our samples are selected at infrared wavelengths and the 
SiR-Fd relation may not stand for galaxies fainter in the in- 
frared or with a low metallicity. However, it is unclear whether 
we should expect higher or lower dust temperatures for these 
types of galaxies because two factors may be in competition. 
On one hand, the calibration of the SFR-to-LiR ratio in Eq. 
assumes an optically thick dust model. The stellar emission is 
strongly attenuated in the galaxies of our samples, as shown by 
the IR to B-band luminosity ratio in Fig.|5]g and this assump- 
tion is justified for our samples. However, galaxies for which 
only a fraction of the intrinsic young stellar emission is atten- 
uated by the dust would have higher SFR-to-LiR ratios (y>l) 
which would give in Eq.|22]lower dust temperatures for a given 
IR surface brightness. On the other hand, for example, Wilson 
et al. ( 1991 1 found that the SFR estimated by the IR luminosity 
was greater than the one estimated by the (extinction-corrected) 
Hq, emission line in the coolest region of the spiral galaxy M 33 
and attributed this excess to the presence of interstellar cirrus 
(Helou [1986:). They can represent a non-negligeable fraction 
(up to 50%-70% or more, Lonsdale-Persson & Helou 1987 
Bell l2003t of the total IR emission in quiescent galaxies. An in- 
crease of the interstellar cirrus contribution would decrease the 
SFR-to-LiR ratio (y< 1 ) that would in turn increase the dust tem- 
perature for a given IR surface brightness. However, it should 
be noted that the most quiescent galaxies in our samples do not 
depart from the £rR-7d relation so it is possible that the two 
factors compensate themselves or have second order effects. 



Table 3. Dependance of the EiR-7d relation on the emis- 
sivity index ft. The effective dust temperature is given by 
7d = a Ej R with Eir being in unit of L pc~ 2 . The empiri- 
cal relation is fitted to the 60/im sample. The model is de- 
scribed in Sect. l6.2l and assumes an absorption mass coefficient 
k(125 yum) = 1.9 m 2 kg -1 and a constant dust-to-gas mass ratio 
?7d = 1/350. 





Observation 


Model 


p 


a b 


a b 


1 


24.32 0.0549 


22.77 0.0571 


1.5 


22.96 0.0503 


22.94 0.0519 


2 


21.68 0.0466 


22.88 0.0476 



It should also be noted that the dust in this model is in ther- 
mal equilibrium and as a consequence, the stochastic heating of 
small dust grains is not taken into account. However, with our 
choice of observables, this fact is mitigated by the adopted def- 
inition of Lir ("Sect. l2~2l which also discards the mid-IR excess 
emission expected from the smaller dust grains. 

7.2. Emissivity index dependence 

The modelled XiR-7d relation (Sect. 16. 2> has a weak depen- 
dence on the emissivity index ft as shown in table[3]by the mod- 
est variation of the scaling factor a and the power-law index b 
given by the relation 7d = a 2^. The modelled values of a and 
b have been calculated with the same values of the absorption 
mass coefficient k{vq) and dust-to-gas ratio rj^ as in Sect. 16.21 
(ft-1.3). Because of the uncertainties associated with these two 
parameters, the comparison between the empirical and mod- 
elled scaling factors a is not insightful, as already noted. More 
noteworthy, the empirical and modelled power-law indices b do 
not depart from each other as ft varies, which makes the agree- 
ment of the model with the observations not sensitive to our 
initial choice of emissivity index. 

7.3. The L IR -r d degeneracy 

With conservative errors of 10 and 15 per cent for the IRAS flux 
densities at 60 and 100/im, the statistical standard deviation of 
the effective dust temperature of the 60 fim dataset is expected 
to be 2.8 K . It is compatible with the sample standard devia- 
tion, which is lower than 3.2 K along the £iR-r d sequence. It 
implies that the introduction of the size parameter has mostly 
disentangled the degeneracy of the LrR-7d diagram down to 
current observational precision. The FIR all-sky survey by the 
satellite Akari will be able to address further this issue. 

7.4. A starburst temperature limit? 

One can wonder what part of the 2^-7^ diagram would be 
populated by ultra-compact starbursts not resolved or too faint 
to be included in our samples and for which the assumption 
of FIR-thin molecular clouds may be incorrect. As the distri- 
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bution of IR emission becomes more compact, the FIR opac- 
ity would increase until the source becomes a blackbody, for 
which the IR surface brightness is related to the temperature by 
Td - (2 IR /4cry5 = o) 1/4 , or numerically by 

/ E \°' 25 

T A = 1.15 ^— K (24) 

Ultra-compact starbursts would lay in Fig. on the left-hand 
side of the steep Eq. |24] and on the upper side of Eq. |23] 
In this scenario, the effective dust temperature would not be 
constrained. However, the FIR-opacity of the most active star- 
bursts is still a matter of debate. Another scenario proposed by 
Klaas et al. (j2001 1 for ultra-luminous galaxies driven by star- 
formation, is that such galaxies are still mostly FIR-thin and 
as a consequence, they would still follow the relation given by 
Eq. |53 until they reach the empirical starburst intensity limit 
by Meurer et al. ( fl997t of 2.0 x 10 5 L pc~ 2 (with a factor of 
3 uncertainty). With these assumptions, the effective temper- 
ature would then be limited to around 44 K for 90% or more 
starbursts. 

8. Conclusions 

We constructed two well-defined local datasets selected at 
60 fim and 170^/m, which are made of 430 IR-bright, non- 
AGN galaxies with reliable radio continuum sizes. We used 
them to investigate the relation between three global parame- 
ters, namely the infrared luminosity L IR , the effective dust tem- 
perature 7^ and the size of the star-forming region r IR . We show 
that 

1 . The dispersion in the Lir-T^ diagram can be explained 
by introducing the size of the star-forming region. 

2. Infrared bright non-AGN galaxies form a plane in the 
(Lir, T&, r IR ) space akin to the fundamental planes of the spiral 
and elliptical galaxies. 

3. The effective dust temperature is related to iW?f R , i.e. 
the IR surface brightness, by a power-law over 5 orders of mag- 
nitude. 

4. Unlike the L lR -T,± relation, the I. lR -Td relation does not 
depend on the IR wavelength used to select or detect galaxy 
samples. 

5. The empirical relation is in agreement with a simplified 
model made of isothermal molecular clouds which are opaque 
in the optical and transparent in the FIR and for which we as- 
sumed a constant emissivity index, gas-to-dust ratio and mass 
absorption coefficient. 

6. The model for which the dust is distributed in a single 
shell around the central starburst is ruled out. 

7. Because 7d also traces the SFR per unit of gas mass, the 
SFR per unit of gas mass correlates with the SFR per unit area. 

8. The infrared compactness turns out to be a parameter 
able to describe the smooth sequence ranging from quiescent to 
starburst galaxies in which the gas surface density, the effective 
dust temperature, the SFR per unit of gas mass and the SFR per 
unit area increase together. 

Assuming that the relation holds for distant galaxies, they 
may be significant, since the hierarchical framework of struc- 
ture formation predicts a decrease of sizes with redshift. It 



would induce an evolution of galaxy colors, which may sta- 
tistically affect the derivation of the cosmic star formation rate 
from infrared galaxy number counts. Such a study will be car- 
ried on a subsequent paper. 

The unprecedented angular far-infrared resolution of ESA's 
HERSCHEL space observatory will allow us to further probe 
the infrared compactness-temperature relation. In a more dis- 
tant future, high redshift galaxies will be spatially resolved by 
missions included in the ESA Cosmic Vision Programme such 
as the Far InfraRed Mission (FIRM) or in the NASA space sci- 
ence roadmap such as: the Single Aperture Far-IR telescope 
(SAFIR), the Space IR Interferometric Telescope (SPIRIT) and 
the Submillimeter Probe of the Evolution of Cosmic Structures 
(SPECS). 
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